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An ordered d-density wave (DDW) state has been proposed as an explanation of the pseudogap 
phase in underdoped high-temperature superconductors. The staggered currents associated with this 
order have signatures which are qualitatively different from those of ordered spins. We apply the 
order parameter theory to an orthorhombic bilayer system and show that the expected magnitude 
as well as the momentum, energy, and polarization dependence of the consequent neutron scattering 
is consistent with the findings of a recent experiment. 

PACS numbers: 71.10.Hf, 75.10.Lp, 75.30.Fv, 71.27.+a 



a. Introduction. The pseudogap phase of the under- 
doped high-temperature superconducting cuprates is no- 
table in its departure from the behavior of a conventional 
metal. Recently, it was proposed that this unusual 
behavior is due to a new broken symmetry, termed d- 
density wave (DDW) ||] , which competes with supercon- 
ductivity. The crucial feature of this order is a staggered 
pattern of orbital currents that spontaneously appear at 
the onset of the pseudogap phase. The underdoped su- 
perconducting state coexists with this order. This pro- 
posal of the existence of an order parameter is funda- 
mentally different from the proposal of staggered current 
fluctuations Q. 

We discuss here the experimental signatures of this 
unusual state vis-a-vis a recent neutron scattering ex- 
periment in underdoped YBC06.6 11 ■ This experiment 
has identified rods of Bragg scattering (within the energy 
resolution window of 1 meV) arising from small magnetic 
moments which increase in strength below the pseudogap 
temperature, with a further increase below the supercon- 
ducting transition temperature. The in-plane wavevector 
is Q = (ir/a, it /a), where a is the in-plane lattice spacing, 
and the intensity is at the level of the background for 
perpendicular wavevector transfer q z = 0, implying an 
"antifcrromagnetic" coupling between the layers within 
a bilayer. Both of these features are precisely the same 
as those observed in the undoped antiferromagnet, where 
the signal is from ordered Cu moments. Such scattering 
from ordered Cu moments deep in the superconducting 
state would be surprising to say the least. 

However, there are four important ways in which these 
experimental results differ from the corresponding ones 
for the undoped antiferromagnet: (1) the magnitude 
of the observed moments is of order 2 x 10 -2 /ib at 
10 K, which is 50 times smaller than that observed in 
the undoped antiferromagnet; (2) the intensity decreases 
rapidly with scattering wavevector q in a manner which 
is inconsistent with ordered Cu spins; (3) aside from the 
elastic peak, there is no significant intensity up to 20 
meV. In other words, there is no evidence for the Gold- 
stone modes (magnons) which must be present if spin- 
rotational symmetry is broken. (4) Finally, it is found 



that there is no 3D order down to the lowest temperatures 
studied despite the fact that the in-plane spatial corre- 
lations are resolution limited, with a correlation length 
greater than 200A. 

We argue that these four are the most robust features 
of the experiment and that they fit the explanation of 
DDW order but are inconsistent with spin order. In or- 
der to understand the experiment, we first introduce the 
theoretical framework of the DDW. Next, we address the 
ratio of the intensities for various Bragg reflections and 
then examine the ratio of the spin-flip to non-spin flip 
scattering for polarized neutrons. 

b. Currents associated with DDW order. The DDW 
order parameter || is a spin singlet particle-hole con- 
densate given by the equal time correlation function 
(a, (3 are spin indices): (c^k+QaCk^) = i$Q Y(k) 6 a/ 3, 
where $q is a real order parameter, and cp is an elec- 
tron destruction operator. The internal degree of free- 
dom of the condensate is given by the d x i_ y i function 
Y(k) = i (cos k x a — cos k y a) . Charge density is not mod- 
ulated in the DDW because Y(k) = 0! It is actually 
the current density that is modulated, as we discuss be- 
low. A conventional charge density wave would occur if 
the internal orbital degree of freedom of the particle-hole 
condensate were of s-wave character, e.g. if Y(k) were 
a constant. The terminology is a very convenient way 
of classifying the internal orbital degree of freedom of a 
particle-hole condensate. 

Neutron scattering from the DDW state is determined 
by its associated current distribution. This distribution 
depends on the details of the current paths along which 
electrons move in a copper-oxide superconductor. As we 
discuss below, unpolarized neutron scattering is fairly in- 
sensitive to these details, but the ratio of non-spin-flip to 
spin-flip polarized neutron scattering is not. 

The charge and spin distribution of an electron at a 
Cu site is given by the well-known Cu form factor (see 
below). However, due to mixing of different orbitals, the 
profile of current flow is more complicated; it will be 
sensitive to many-body effects when the current is sub- 
stantial, and it will be affected by the bi-layer coupling 
and, in YBCO, by the influence of the chains. We make 
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a simple and tractable Ansatz for this current distribu- 
tion, including the effects of orthorhombicity due to the 
chains. The orthorhombicity, even though it is small in 
terms of the lattice constants, is well known to give rise to 
remarkable anisotropy; for example the a — 6-anisotropy 
of the superfluid density can be as large as 2-3 || , and 
perhaps even larger 0. 

Thus, we define the Fourier transform of the expecta- 
tion value of the current, (j(q)), satisfying q- (j(q)) = 0, 
by 
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(2TrH/a,2TrK/a 7 2irL/c) is a magnetic reciprocal lattice 
vector, where the lattice constants are a = 3.86A and 
c = 11.82A; we ignore slight orthorhombicity in the 
lattice constants, but not in the physical quantities in 
which it is magnified. Note that /(q) = sin(2|^), where 
d = 3.25A is the intrabilayer separation. The temper- 
ature and doping dependence of the currents are deter- 
mined by the DDW order parameter $q, which must 
appear as a thermodynamic phase transition. 

If the electrons behaved classically and took straight 
line trajectories between point-like Cu sites, we would 
have a(q) = /3(q) = 1. In reality, however, these are 
non-trivial functions of q. Guided by the phenomenol- 
ogy, we choose the simplest model in which these are 
non-trivial functions of q z alone. At long scales, the cur- 
rent distribution is essentially classical, a(0) = (3(0) = 1. 
At short distances, the thickness of the current trajecto- 
ries becomes apparent. We assume that currents flowing 
along the y axis are more spread out, reflecting the un- 
derlying orthorhombicity of YBCO, a(q z ) 7^ (3(q z ) for q z 
large. In tetragonal materials, the asymmetry could be 
due to spontaneous symmetry breaking |1. 

The requirement a(0) = (3(0) — 1 automatically leads 
to zero average current along z in the bilayer model be- 
cause of the bilayer form factor /(q). If we want to inter- 
pret a(q z ) and (3(q z ) to be merely form factors and want 
to adhere to the original order parameter we can 

do so if we stipulate that a(q z ) — (3(q z ) ~ q z , as q z — > 0. 
This will make the average z-component of the current 
vanish (as it must) for the single layer model for which 
the same expression for the current holds, except that 
the factor /(q) is missing. 

The orbital currents give rise to magnetic fields of order 
10 G, and a very small magnetic moment, of order 

4 x 10 hb, consistent with the experimental estimate of 
order 2 x IO^b ©■ 

c. Unpolarized elastic neutron scattering. We first re- 
call spin-scattering and then discuss scattering from or- 
bital currents. The well-known scattering cross section of 



unpolarized neutrons from localized Cu spins in YBCO, 
when the net spin moment is along the c-axis, is 
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where q is the momentum transfer, and the atomic form 
factor g(q) is the Fourier transform of the normalized 
density of unpaired electrons in a single ion pj. We 
have assumed that the spins are antiferromagnetically 
coupled within an YBCO bilayer and that different bi- 
layers are uncoupled, so that the scattering cross-section 
is proportional to that due to a single bi-layer. 

Similarly, there is a formula which applies when the 
spin moment lies in the x — y plane ]lC[ |. If the spin 
direction in the plane is averaged over magnetic domains 
or averaged due to spin glass order, this is 
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It is worth emphasizing that the quantity in large paren- 
theses in (|^) increases with q z , as does its counterpart in 
the case in which the spins point in some direction in the 
a — b plane which is not averaged over. 

Next, we consider Bragg scattering from orbital cur- 
rents M. The neutron magnetic moment generates a 

(r — r ) 

vector potential, A, given by A = /j x | r °-r % > wnere 
fi = — 1.91(^-)s„, s„ is the neutron spin, and m„ is the 
mass of a neutron; the electron and the neutron coordi- 
nates are r e and r ra , respectively. The coupling of the 
electrons to this gauge field is given in momentum space 
by the Hamiltonian H int = J j(q) • A(q). The Bragg 
intensity of unpolarized neutrons is then 
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where X(q z ) — a(q z ) / ' (3(q z ) . We emphasize that there 
are no approximations in this formula, such as the SU(n) 
mean field approximation ^| ; it simply follows from the 
assumption of DDW order. 

We estimate the factor (3 2 (q z ) by g 2 (q) for Cu spins 
found by Shamoto et al. fic|| as an upper bound, because 
(3 2 (q z ) must fall off faster than g 2 (q) as the charge distri- 
bution for orbital current loops is more spread out than 
the atomic orbitals. The case for orbital currents made 
below will be even stronger with the true form factor. 

The most salient feature distinguishing (||) and (||) 
from (ji)) is apparent from dimensional analysis. The for- 
mer scale with q as q° while the latter dominantly scale 
as l/g 4 (recall that |j(q)| ~ l/q according to (|j) and 
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a{q z ) and /3(q z ) are of order unity). The latter form fol- 
lows from the inherent spatial extent of a current loop 
and leads to a rapid decay of the scattering cross-section 
with q. 

Let us consider the ratios of the Bragg scattering in- 
tensities at various magnetic reciprocal lattice vectors; 
we will quote the intensities at the representative values, 
X(q z ) = 1 and 3, thereby crudely replacing X(q z ) by a con- 
stant for wavevectors larger than or of order 2ir/c. From 
(|]), we see that the ratio of the intensity at (0.5, 0.5, 1) 
(in reciprocal lattice units) to that at (1.5,1.5,1) is 111 
at A = 1, and 23 at A = 3. Meanwhile in the spin case, 
the same ratio is 1.9 if the spins lie in the a — b plane and 
1.4 if they are along the c-axis. The strong suppression 
of the intensity at (1.5, 1.5, 1) in the orbital case is con- 
sistent with experiments [p|]l l|] where no peak is observed 
at (1.5, 1.5, 1) or (1.5, 1.5, 2), while peaks are observed at 
(0.5, 0.5, 1) and (0.5, 0.5, 2). The absence of an observable 
peak at (1.5, 1.5, 1) is inconsistent with spin order. 

Let us now contrast the ratios of the intensities at dif- 
ferent values of q z for spins pointing within the plane, 
Eq. |[ and orbital currents, Eq. ||, with the form factor 
correction described above. The ratio of the intensity at 
the Bragg reflection (0.5, 0.5, 1) to that at (0.5, 0.5, 2) is 
0.96 for orbital moments at A = 1, increasing to 1.5 at 
A = 3. For spins pointing in the plane it is 0.5. Experi- 
mentally H it is found to be 1.16±0.16, suggesting orbital 
currents as the source of scattering. One can easily see 
from the known Cu form factor Jl0| that the correction 
due to the atomic form factor for nearby Bragg reflec- 
tions, such as (0.5,0.5, 1) and (0.5,0.5,2), is negligible. 

Moreover, for spin-moments lying parallel to the plane, 
the ratio of the intensity at (0.5,0.5,5) to (0.5,0.5,1) 
should be 1.6. This is clearly not the case in the ex- 
periment [§,0 because there is no intensity above the 
background at (0.5, 0.5, 5). For orbital currents, the same 
ratio is 0.2 (which is only an upper bound) at A = 1 and 
0.1 at A = 3; as a result, there may not be a measur- 
able intensity at (0.5, 0.5, 5) as discussed in Ref. ||, thus 
confirming once again the source of the scattering to be 
orbital currents. 

d. Polarized elastic scattering. For polarized neu- 
trons, we have 
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where a) and \(3) are the initial and final spin states of 
the neutrons. It is easy to see that if the neutron polar- 
ization is parallel to the scattering vector, the scattering 
is entirely spin flip irrespective of (j(q)). 

In contrast, if it is perpendicular to the scattering vec- 
tor, then the scattering can be either spin flip or non-spin 
flip depending on (j(q)} and q. If q = (H, H, L) and the 
neutrons are polarized perpendicular to q, in the [1, 1,0] 
direction, the ratio of the non-spin-flip to spin flip scat- 
tering intensities is: 



While this ratio vanishes at A = 1, it increases rapidly 
with A because c ^> a. At (0.5,0.5,1) and A = 3 it is of 
order unity, more precisely 1.4, which is similar to the 
ratios found in the experiments of Refs. |0,Q. 

e. Interlayer Coupling. Earlier, we assumed that the 
order parameters from the two layers in a bilayer are 
opposite in sign. Let us see why this should be so. In- 
teractions which couple the density at a site in one layer 
to the density at a site in another layer will not cou- 
ple the DDW order parameters in the two layers because 
they will average over directions in the a — b plane, giving 
zero. However, interlayer tunneling will couple the DDW 
order parameters at second order. 

Consider the interlayer tunneling Hamiltonian []l3| be- 

tween layers (1) and (2): H c = - £ kct fi_(k)(c£2 f c£> + 

h. c), where t±(k) — (tj_/4)(cos k x a — cosk y a) 2 . Let \s) 
be the ground state for t± = with equal DDW order 
parameters in each layer, while |o) is the ground state 
for tj_ = with DDW order parameters equal and op- 
posite. The energy difference between these two states, 
per site, is given to lowest non-vanishing order in t± by 

V » -l£ $ Q- where A n = V(A DDW ) 2 + (A DSC ) 2 , 
and A DSC is the maximum of the superconducting gap, 

i. e. the antisymmetric state has lower energy. The above 
result also holds for body-centered tetragonal materials, 
with the modified tunneling matrix element [Q. More- 
over, we see why the coupling between bilayers is so 
small: the tunneling matrix elements between bilayers 
are much smaller than those within a bilayer because 
c = 11.82A > d = 3.25A. 

Because the coupling between the unit cells is excep- 
tionally small, we expect to see 2D Bragg rods without 
ZD order for temperatures which are not exceedingly low; 
recall that the transition to the DDW state is an Ising 
transition unlike the corresponding spin problem. In- 
deed, as remarked earlier, Mook et al. [[| find such Bragg 
rods. 

/. Inelastic scattering. We note that because DDW 
order is Ising-like, there are no Goldstone modes, but a 
gap in the frequency spectrum. This is consistent with 
the experimental finding of a spin gap above the elastic 
peak H , extending up to an energy of 20 meV. For spin- 
ordering, a Bragg peak will necessarily lead to a Gold- 
stone mode (the Ising anisotropy in these materials is 
negligibly small; see, for example, Ref. po[.), which is 
continuously connected to it. It is not possible to have 
one without the other. 

The magnetic resonance peak [ p"5| in inelastic scatter- 
ing at 34 meV at the same momentum, Q, also testifies 
against the spin-ordering scenario. If the spins are or- 
dered, this resonance should be at zero energy - i.e. it 
should not be distinct from the Bragg peak with the same 
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quantum numbers! Furthermore, in the presence of bro- 
ken spin-rotational symmetry, excitations are not orga- 
nized into SU(2) multiplets (but, rather, U(l) multiplets 
if rotations about one axis are preserved), so the reso- 
nance could not possibly be a triplet mode as suggested 
in Rcf. @. 

g. Global Phase Diagram. According to the interpre- 
tation proposed here, DDW order has been observed in 
neutron scattering in YBC06.6 |§- On the other hand, 
neutron scattering in YBCO6.35 produces a signal which 
is 10 times more intense, is due to moments which lie in 
the a — b plane and is believed to be due to spin glass 
ordering ||. This not only supports the hypothesis that 
the effect observed in YBCOg.6 is due to something other 
than spin ordering - namely DDW ordering - but also 
suggests that DDW order disappears at low doping in 
favor of spin glass order and, possibly, other competing 
orders as well || . 

In neutron scattering experiments on YBCO6.5, Sidis 
et al. [fl2f find similar magnetic Bragg peaks with a ratio 
of 0.67 between the intensity at (0.5, 0.5, 1) and that at 
(0.5,0.5,2). This is intermediate between the observed 
ratio in the putative DDW state of YBCO6.6 and the 
observed ratio in the spin glass of YBCO6.35. The mea- 
surements of Ref. [[l2| seem to be indicative of the spin 
glass phase. It is clear therefore that the doping depen- 
dence may be crucial, changing the results as we move 
from YBC0 6 . 5 to YBC0 6 . 6 . 

With the above ideas as inspiration, we propose the 
global phase diagram of Fig. [l| 
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FIG. 1. Phase diagram. 3D AF is the three-dimensional 
antiferromagnetic phase, and TV is the corresponding Neel 
temperature. Competing phases refer to a complex set of 
charge-ordered states; tcp is a tetracritical point; x c is a quan- 
tum critical point, and T* signifies the pseudogap transition. 

h. Summary. To summarize, the ratio of spin-flip to 
non-spin-flip scattering depends on details of the current 
distribution, but the rapid decrease of the scattering in- 
tensity with q is a robust feature which follows from 
the inherent spatial extent of the current loops (which 
is much larger than a Cu or even an O orbital). We 
believe that DDW order is the only way of reconciling 



the constraints following from the wavevector and polar- 
ization dependence of the neutron scattering intensity. 
These phenomenological considerations, although quite 
tight, clearly await a microscopic understanding. 

Note added: Recent striking zero-field muon spin re- 
laxation measurements in YBCOg.67 and YBCO6.95 have 
uncovered the onset of static magnetism consistent with 
our DDW picture, and the phase diagram proposed above 
]l6[ , in particular, a phase transition within the super- 
conducting dome. 

S. C. is supported by NSF-DMR-9971138, and C. N. 
by NSF-DMR-9983544 and the A.P. Sloan Foundation. 
The work of H. -Y. K. was conducted under the aus- 
pices of the DOE, supported by funds provided by the 
University of California for the conduct of discretionary 
research by Los Alamos National Laboratory. We thank 
E. Abrahams, B. Keimer, S. Kivelson, R. B. Laughlin, 
and J. Tranquada for discussions. We especially thank 
H. Mook and J. E. Sonier for interesting and important 
correspondence. 



[1] 

[2] 

[3] 

[4] 
[5] 

[6] 
[7] 
[8] 



[9] 

[10] 
[11] 

[12] 
[13] 

[14] 
[15] 



[16] 



S. Chakravarty, R. B. Laughlin, D. K. Morr, and C. 

Nayak, Phys. Rev. B 64, 094503 (2001). 

S. Chakravarty and H.-Y. Kee, Phys. Rev. B 61, 14821 

(2000) . 

C. Nayak, Phys. Rev. B 62, 48 80 (2000); C. Nayak and 
F. Wilczek, |cond-mat/9510132i 

P. A. Lee arid X. -G. Wen, |cond-mat/00084li . 

H. Mook, P. Dai, and F. Dogan, Phys. Rev. B 64, 012502 

(2001) . 

D. Basov et al, Phys. Rev. Lett. 74, 598 (1995) 
A. G. Sun et al, Phys. Rev. B 52, R15731(1995). 

S. Kivelson, E. Fradkin and V. J. Emery, Nature 393, 
550 (1998); see also, S. Caprara, et al., J. de Phys. IV 
(Proceedings) 9, 329 (1999); J. Zaanen, Physica C 317- 
318, 217 (1999); J. M. Tranquada, et al, Nature 375, 
561 (1995), and Phys. Rev. B 54, 7489 (1996). 
T. C. Hsu, J. B. Marston, I. Affleck, Phys. Rev. B 43, 
2866 (1991). 

S. Shamoto et al, Phys. Rev. B 48, 13817 (1993). 
H. A. Mook, Pengcheng Dai, S. M. Hayden, S. -H. Lee, 
A. Hiess, and F. Dogan, preprint, submitted to Science. 
Y. Sidis et al, Phys. Rev. Lett. 86, 4100 (2001). 
S. Chakravarty, A. Sudbo, P. W. Anderson, and S. 
Strong, Science 2 61, 337 (1993). 

E. Pavarini et al. 



:ond-mat/0012051 



See, for example, P. Dai et al, Phys. Rev. B 63, 054525 
(2001) and H. F. Fong et al, Phys. Rev. B 61, 14773 
(2000). 

J. E. Sonier et al, Science 292, 1692 (2001). 



4 



